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PREFACE  •  This  work  was  undertaken  in  connection  with  studies  to 
develop  theories  and  computational  methods  for  scattering  of  radar  waves 
hy  vegetation.  To  accomplish  this,  radiative  transfer  theory  has  been  extended 
to  layered  media,  where  the  scattering  parameters  vary  from  layer  to  layer 
and  where  the  scattering  parameters  for  any  given  layer  can  be  computed  from 
a  knowledge  of  the  distributions  of  vegetative  scattering  elements  and  ap¬ 
propriate  theories  for  single  scattering. 

I  he  extended  theory  eventually  became  a  multi-channel  (three-dimensional) 
formulation.  The  first  steps  in  deriving  this  theory  followed  the  historical 
precedent  ot  deriving  the  formulas  to  treat  the  one-dimensional  case,  which 
is  treated  in  the  present  report.  Once  this  pattern  was  set.  the  equations  were 
further  extended  to  cover  the  multi-channel  case,  and  these  results  will  be 
given  in  a  separate  report. 


In  addition  to  setting  a  mathematical  pattern  to  be  followed  in  subsequent 
work,  a  case  tor  diffuse  reflection  was  developed  from  a  semi-infinite  homog¬ 
eneous  medium  by  using  a  mathematical  limiting  process.  Since  the  results 
are  ot  general  interest,  reflection  coefficients  were  computed  and  tabulated 
tor  a  wide  range  ot  values  ot  the  albedo  for  single  scattering. 


This  work  was  performed  in  the  Center  for  Theoretical  and  Applied  Physical 
Sciences  ot  the  Research  Institute  ot  the  U.S.  Army  Engineer  Topographic 
Laboratories  under  work  unit  4A161 I02B52C  S3  0003. 


C  OL  Daniel  L .  Lycan.  Ch  was  Commander  and  Director  and  Mr.  Robert  P. 
Macchia  was  Technical  Director  of  ETL  during  the  report  preparation. 
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RADIATIVE  TRANSFER  IN  ONE-DIMENSIONAL, 
DISCRETELY  STRATIFIED  MEDIA 


INTRODUCTION  •  I  his  research  report  discusses  one-dimensional 
transfer  of  radiation  in  a  medium  composed  of  layers  whose  scattering  pro¬ 
perties  vary  from  one  level  to  the  next.  In  addition,  it  is  assumed  that  no 
radiation  sources  are  present  within  the  medium.  The  case  considered  is  for 
diffuse  reflection  and  transmission  of  externally  incident  radiation.  The  scat¬ 
tering  properties  of  a  single  layer  are  assumed  to  he  known,  and  the  math¬ 
ematical  relations  are  developed  for  expressing  the  scattering  properties  of  the 
multi-layer  in  terms  of  the  scattering  parameters  of  the  single  component 
layers. 


This  work  was  undertaken  to  determine  radar  scattering  by  vegetation.  Since 
rigorous  mathematical  solution,  are  generally  known  only  for  homogeneous 
regions  of  relatively  simple  goemetry.  it  was  considered  appropriate  to  use 
a  discrete  layer  approach  to  account  for  the  variation  of  scattering  parameters 
with  depth  inside  the  vegetation.  The  radiation  fields  have  been  taken  to  be 
one-dimensional  to  establish  the  pattern  of  the  mathematical  solution.  Sub¬ 
sequently,  the  solutions  have  been  extended  to  radiation  fields  of  higher 
dimensionality  involving  a  discrete  channel  model.  The  results  of  the  solutions 
will  be  given  in  separate  reports,  with  the  present  report  restricted  to  the 
simplified,  one-dimensional  case.  This  method  of  development  is  consistent 
with  the  historical  precedent  set  by  Schuster.1  whose  solution  of  a  one-dimen¬ 
sional  problem  is  generally  considered  to  have  lain  the  foundation  for  all  fur¬ 
ther  developments  of  radiative  transfer  theory. 

In  addition  to  establishing  a  framework  upon  which  to  build  the  discrete  multi¬ 
channel  theory,  this  research  report  also  gives  a  new  one-dimensional  solution 
lor  back  reflection  from  a  homogeneous  semi-infinite  medium.  To  obtain 
this  solution,  an  invariance  principle  is  applied  to  the  transfer  equations  in 
their  discrete  lorm.  I  lie  transition  to  the  case  of  a  continuous  scattering 
medium  is  then  obtained  by  a  mathematical  limiting  process.  Since  the  re¬ 
sulting  equation  lias  a  simple  lorm.  a  table  of  backscattering  functions  has 
been  computed  and  tabulated. 


Schuster.  Arthur.  Radiatum  Thrmnih  u  /nwi'  \minspfnrt  ,  ||»c  Astiopln sU.il  Joiirn.il.  Vol  2\. 
No.  I.  Kimiurv ,  1905,  pp.  I -22. 
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GENERAL  FORMULATION  •  The  scattering  properties  of  a  single 
component  layer  will  be  specified  by  its  scattering  coefficient  a  and  its 
transmission  coefficient  r.  In  practice,  these  coefficients  may  be  taken  from 

a.  Solutions  of  multiple  scattering  theories  of  radiative  transfer  theory 

b.  Solutions  of  single  scattering  theories,  if  the  layer  is  thin  enough 
that  the  amount  of  multiple  scattering  within  the  layer  is  negligible 

c.  Composite  layers  that  have  been  assembled  from  layers  of  type 
a  or  b,  using  the  methods  to  be  presented  below 

The  intensities  of  radiation  flux  in  the  forward  and  backward  directions  will 
be  specified  by  I  and  R.  respectively  (figure  1). 


i 

t 


FIGURE  1.  Single  Component  Layer  Scattering  Parameters. 
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In  radiative  transfer  theory',  it  is  customary  to  use  the  term  “scattering 
(or  reflection)  function  (or  coefficient)” to  specify  a  ratio  of  a  backscattered 
intensity  to  an  incident  intensity.  Similarly,  the  term  “transmission  function 
(or  coefficient)” is  used  to  denote  the  ratio  of  a  transmitted  intensity  to  an 
incident  intensity.  These  coefficients  will  be  defined  more  precisely  below. 
In  general  it  is  necessary  to  use  different  coefficients  for  the  scattering  and 
transmission  functions  in  opposite  dilections.  However,  for  thin  single-scat¬ 
tering  layers,  symmetry  properties  of  the  scatterers  and  reciprocity  theorems 
for  single  scattering  can  sometimes  be  used  to  demonstrate  that  the  same 
coefficients  can  be  used  for  both  forward  and  backward  directions  under 
certain  conditions.  F.ven  when  the  reciprocity  holds  for  single  layers,  it  may 
fail  for  thicker  multi-layers.  This  phenomemon  is  known  as  polarity.  An 
example  of  its  occurrence  would  be  a  double  layer  with  one  half  purely  ab¬ 
sorbing  and  the  other  half  purely  scattering. 


For  the  single  component  scattering  layer  being  considered,  the  incident 
fluxes  are  I  and  R,  .  The  subscripts  <>  and  I  denote  that  the  radiation 
fields  lie  respectively  to  the  left  or  right  of  the  layer.  In  general,  the  scat¬ 
tering  and  transmission  coefficients  depend  on  whether  the  incident  intensities 
are  in  the  forward  or  backward  direction.  Thus,  if  R,  =  0.  then  a  and 
r  are  defined  by  a  =  R  /Io  and  t  =  I(  /l().  On  the  other  hand,  if 
l()  =  0,  then  o'  =  l(.  /R(  and  t  =  Ro /R,  .  In  the  general  case,  neither 
lo  nor  R,  are  assumed  to  be  zero.  Then,  the  emergent  fields  !,  and 
R()  are  given  by  the  linear  superposition  of  the  contributions  from  both 
incident  fields  1  and  R,  . 

O  I 


1,  =  rl  +  o’R 

1  o  I 


( 1 1 


R  =  ol  +  r'R. 

o  o  I 


(?) 


In  the  case  of  thin  single-s  altering  layers,  reciprocity  theorems  that  apply 
to  single  scattering  may  be  equivalent  to  the  condition  a  =  o'  and  t  =  r'. 
It  will  be  shown  in  the  following  section  that  when  two  different  layers  are 
combined,  the  corresponding  coefficients  for  the  composite  layer  will  generally 
no  longer  satisfy  such  a  reciprocity  condition. 


COMBINATION  OF  LAYERS  •  The  relationship  lor  scattering  and 
transmission  by  two  adjacent  layers,  whose  individual  scattering  and  trans¬ 
mission  coefficients  are  given,  will  now  be  derived.  The  geometrical  arrange¬ 
ment  is  shown  in  figure  2. 


FIGURE  2.  Combination  of  Two  Component-Scattering  Layers. 


The  intensities  of  the  radiation  fields  between  the  layers.  IM  and  RM . 
will  contain  not  only  transmitted  radiation  but  also  scattered  radiation.  There¬ 
fore.  the  possibility  of  multiple  scattering  between  layers  now  arises,  which 
will  be  reflected  in  the  resulting  equations. 

If  the  left  layer  is  considered  by  itself,  the  emergent  intensities.  IM  and  R(>. 
may  be  expressed  in  terms  of  the  incident  intensities.  l(i  and  RM  .  in  the 
manner  given  by  equations  1  and  2. 
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(4) 
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When  equations  9  and  10  are  used  to  eliminate  l%1  and  RM  from  equation 
3  after  some  algebraic  manipulations,  it  is  found  that  l(  can  be  expressed 
in  terms  of  l(). 


I, 


tr  7. 


“  °i  ,jr 


(11) 


Also,  when  equations  9  and  1  1  are  combined  with  equation  4  to  eliminate 
RS1  and  1,  ,  it  is  possible  to  express  R  (  in  terms  of  I 


By  an  analogous  procedure,  if  10  =  0.  so  that  the  only  source  present  is 
Rj. .  equations  3  and  4  become 


which  can  be  solved  for  RM 


and  I.,  in  terms  of  R,. 

M  t> 


Rv,  =  rj  1  R 

M  L  o 


'm  = 


(15) 

(16) 


By  inserting  equations  15  and  16  into  equation  (i  to  eliminate  RM  and 
IM  ,  an  expression  is  found  for  R(>  in  terms  of  R,  . 


When  equations  16  and  17  are  combined  with  equation  5  to  eliminate 
IM  and  Ro  ,  an  expression  for  Ij  in  terms  of  R,,  is  found. 


(18) 


The  total  solution,  when  i0  and  R,  are  both  non-zero,  is  given  by  the 
superposition  of  the  solutions  given  in  equations  11,  12,  17  and  18. 


(19) 


(20) 


liquations  19  and  20  give  the  emergent  fields.  1,  and  R  ,  for  the  double 
layer  in  terms  of  the  incident  fields  lo  and  R,  .  By  analogy,  with  the  single 
layer  case  given  by  equations  1  and  2.  it  is  possible  to  characterize  the 
scattering  properties  of  the  composite  double  layer  by  reflection  and  trans¬ 
mission  coefficients  a(. ,  a'c.  r(. .  r'(. . 


=  r  1  + 

C  o 

ac 

R, 

(21) 

R 

O 

=  a  1  + 

(  o 

t 

Tc 

R, 

(22) 

T 

T  T 

R  L 

\ 

1  -  o’L 

aR 

(23) 

i 

T 

_  T'l  rR 

(24) 

V 

1  -  °R 

°L 

10 
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a,  +  'L  --S-Ik 

1  -  o' 


L  °R 


(25) 


a 


c 


(26) 


Equations  23  through  26  show  that  unless  the  single  component  layers 
are  identical  and  have  reciprocal  scattering  properties  wirh  respect  to  direction, 
reciprocity  does  not  hold  for  the  composite  double  layer. 


Clearly  the  formulas  that  have  been  developed  in  this  section  may  be  applied 
repetitively  to  build  a  scattering  model  for  a  medium  consisting  of  many 
layers,  if  the  scattering  and  transmission  properties  of  each  individual  layer 
are  known. 


INTERPRETATION  OF  THE  MULTIPLE  SCATTERING  TERMS  • 

Equations  23  through  26  contain  fractional  terms  that  account  for  the 
presence  of  multiple  scattering.  If  these  terms  are  expanded  in  a  geometric 
series,  each  resulting  term  is  found  to  represent  a  set  of  reflections  and  trans¬ 
missions  by/through  the  single  component  layers.  This  will  be  illustrated  by 
taking  equation  25  as  an  example.  Using  a  geometric  series  transforms  equa¬ 
tion  25  to  the  following  form: 


°c  = 


+  °R  Tl.  +  TL  °R  °L  aR  TL 


+  TL  °R  °l  °R  °L  aR  Tl.  + 


(27) 


The  first  term  accounts  for  radiation  scattered  once  by  the  left  layer.  The 
second  term  accounts  for  radiation  transmitted  through  the  left  layer,  re¬ 
flected  by  the  right  layer,  and  transmitted  back  through  the  left  layer.  In  this 
manner,  the  factors  comprising  any  given  term  account  for  the  history  of 
reflections  and  transmissions  represented  by  each  term.  This  is  demonstrated 
schematically  in  figure  3. 


II 


THE  SEMI-INFINITE  HOMOGENEOUS  CASE  TREATED  BY  AN  INVARI¬ 
ANCE  PRINCIPLE  •  The  reflection  coefficient  for  a  semi-infinite  medium 
composed  of  layers,  all  having  the  same  scattering  characteristics,  can  be 
calculated  by  using  an  invariance  principle.  The  medium  will  be  described 
by  a  notation  similar  to  the  one  used  in  the  preceding  sections,  but  with  a 
new  set  of  subscripts  as  shown  in  figure  4. 


ETC. 


•  #  « 


FIGURE  4.  Reflection  from  a  Semi-Infinite  Stack  of  Scattering  Layers. 


The  reflection  coefficients  for  the  complete  array  ol  scattering  layers  is  given 

by 


R  =  SI  (28) 


and  the  scattering  and  transmission  coefficients,  a  and  t.  are  the  same  for 
each  layer.  Since  the  scattering  array  remains  unchanged  if  the  left  layer  is 
removed,  the  condition  of  invariance  can  be  stated  as  follows: 


R 


l 


SI, 


<2‘» 
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The  emergent  intensities  tor  the  Virst  plate,  expressed  m  terms  ol  the  incident 
intensities,  take  the  form 


R«  =  rK,  +  «'o 


'l  =  Tl„  +  °RI 


and  the  quantities  R0  and  R]  can  he  eliminated  by  means  of  equations 
28  and  2l>. 


S1»  =  rS1,  +  °'0 


=  rl0  +  <*l, 


Both  l0  and  I,  can  be  removed  from  equations  32  and  33.  resulting  in 
a  quadratic  equation  for  S  in  terms  of  a  and  r. 


a S2  +  (t2  -  o2  -  1)  S  +  a  =  0 


Equation  34  has  a  solution  of  the  form 


+  a2  -  t2 


and  equality  holds  only  il  no  absorption  occurs.  Squaring  both  sides  ol  3(> 
(remembering  that  they  are  both  positive)  leads  to  the  condition 


which  combined  with  the  fact  that 


S  <  1 


(38) 


indicates  that  the  negative  sign  must  be  chosen  in  equation  35. 


S  =  X  -  yj X2  -  I 


X 


1 


-) 

+  CT 


1(3 


r“ 


(39) 


(40) 


It  is  of  interest  to  note  t hat  if  equation  34  is  written  as 


t~  S 

S  =  o  +  -  (41) 

1  -  oS 

and  the  second  term  on  the  right  is  developed  as  a  geometric  series,  an  inter¬ 
pretation  similar  to  that  given  in  the  preceding  section  can  be  made.  Also, 
equation  41  can  be  used  to  develop  an  iteration  scheme  in  which  an  approxi¬ 
mate  value  lor  S.  inserted  on  the  right  side,  will  yield  a  better  approximation 
for  S.  These  ideas  take  on  more  importance  in  a  multi-channel  formulation 
and  will  not  be  developed  further  in  this  report. 


PASSAGE  TO  THE  CONTINUUM  AS  A  MATHEMATICAL  LIMITING 
PROCESS  •  The  equations  that  have  been  developed  up  to  this  point 
are  applicable  to  layers  of  finite  thickness  whose  scattering  properties  are 
expressed  by  a  and  t.  The  single-scattering  characteristics  of  the  medium 
must  be  expressed  ideally  in  terms  of  scattering  parameters  that  pertain  to  an 
infinitesimal  layer  as  its  thickness  approaches  zero.  Since  in  this  limit  both 
a  and  t  always  vanish,  they  cannot  be  used.  Since  it  is  customary  to  use  the 
albedo  for  single  scattering  <cr).  an  analogous  quantity  will  be  defined  for  this 
one-dimensional  case. 


In  conformity  with  Schuster.2  it  will  be  assumed  that  half  of  the  scattered 
radiation  will  appear  in  the  backward  direction  and  half  will  appear  in  the 
transmitted  flux.  The  fraction  of  radiation  scattered  is  given  by  the  scattering 
coefficient 


s 


2a 


(42) 


and  the  fraction  removed  from  the  incident  flux  by  both  scattering  and  ab¬ 
sorption  is  given  by  the  extinction  coefficient 


e  =  1  -  (r  -  o) 


(43) 


where  the  term  is  parentheses  accounts  for  the  fact  that  the  transmitted  |lux 
contains  half  of  the  scattered  radiation.  If  the  albedo  for  single  scattering  is 
defined  by 


s 

vr  =  — 
>’  e 


(44) 


it  is  then  possible  to  express  a  and  t  as  follows: 

'ey  e 

a  =  (45) 


1 


"Schuster.  Arthur.  Radiation  Through  a  h'ogj!}'  Atmosphere.  The  Astrophysieal  Journal.  Vol.  21, 
No.  I.  January,  1905.  nn.  1-22. 
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inserting  these  values  into  equation  40  yields: 


2  -  'CJ  -  e  (1  -  “G7  ) 

X  =  - a - si- 

'TD 


Furthermore,  letting  the  layer  thickness  approach  zero  is  mathematically 
equivalent  to  allowing  the  extinction  coefficient  to  approach  zero. 


X  =  lint  X 


Insertion  of  this  value  into  equation  39  results  in  a  simple  expression  for 
the  reflection  coefficient  S  of  a  semi-infinite  homogeneous  medium  whose 
scattering  properties  are  specified  entirely  by 

~  2  -  "ccr  -  2  -J  1  -  txr 


S  =  - 


(I  -  J  1  -  -c 7  1- 


In  practice.  ao  would  be  obtained  from  the  values  of  s  and  e  for  a  thin 
layer  where  the  amount  of  multiple  scattering  would  be  very  small. 


TABULATED  BACKSCATTERING  FUNCTIONS  FOR  A  HOMOGENEOUS 
SEMI-INFINITE  MEDIUM  •  Reflection  coefficients  S  have  been  com¬ 
puted  for  a  wide  range  of  values  of  the  albedo  for  single  scattering  oo  by 
using  equation  50.  Since  they  may  be  of  general  interest,  they  are  given 
below  in  table  1 . 


TABLE  1. 

REFLECTION  COEFFICIENTS  S  (rzr  ) 

O 


0.0001 

0.001 

0.01 

0.05 

0.1 

o.: 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

0.95 

0.99 

0.999 

0.9999 

0.99999 

1.0 


0.000025 
0.00025 
0.00251 
0.01282 
0.02633 
0.05573 
0.08893 
0.1270 
0.1716 
0.2251 
0.2922 
0.3820 
0.5195 
0.6345 
0.8182 
0.9387 
0.9802 
0.9937 
1 .0000 
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CONCLUSION  •  In  order  to  set  the  mathematical  pattern  of  development 
of  a  suitable  theory  of  radiative  transfer  for  use  in  studies  of  radar  scattering 
by  vegetation,  a  theory  of  one-dimensional  radiative  transfer  in  discrete  layered 
media  was  developed.  By  using  this  theory,  it  was  shown  how  single  layers 
with  different  scattering  properties  can  be  combined  to  obtain  the  scattering 
parameters  for  multi-layered  media.  An  interpretation  of  the  resulting  equa¬ 
tions  was  given  that  accounts  for  all  possible  multiple  scattering  between  the 
component  layers. 


A  semi-infinite  homogeneous  medium  was  also  treated  by  using  an  invariance 
principle.  In  this  case,  the  transition  was  made  from  the  discrete  to  the  con¬ 
tinuous  case,  and  a  simple  formula  was  given  for  the  reflection  coefficient 
in  terms  of  the  albedo  for  single  scattering.  Since  the  results  are  of  general 
interest,  a  table  covering  a  wide  range  of  circumstances  was  prepared. 
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